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Core-collapse supernovae (CC-SNe) are the explosions that announce 
the death of massive stars. Some CC-SNe are linked to long-duration 
gamma-ray bursts (GRBs) and are highly aspherical. One impor- 
tant question is to what extent asphericity is common to all CC- 
SNe. Here we present late-time spectra for a number of CC-SNe 
from stripped-envelope steirs, and use them to explore any cispheric- 
ity generated in the inner pcirt of the exploding stcir, necir the site 
of collapse. A range of oxygen emission-line profiles is observed, in- 
cluding a high incidence of double-peaked profiles, a distinct signa- 
ture of an cispherical explosion. Our results suggest that all CC-SNe 
from stripped-envelope steirs Eire aspherical explosions and that SNe 
accompanied by GRBs exhibit the highest degree of asphericity. 

Massive stars (M > lOM©) end their lives when the nuclear fuel in their innermost 
region is consumed; lacking sufficient internal pressure support, they can no longer with- 
stand the pull of gravity. Their core then collapses to a neutron star or a black hole. 
Gravitational energy from the collapse produces an explosion that expels the rest of the 
star in what is observed as a supernova (SN). 

Core-collapse SNe (CC-SNe) are classified (i) by how much of the stellar envelope 
is present at the time of the explosion {2). Stars that retained their hydrogen envelope 
produce SNe with a H-rich spectrum, classified as Type II. On the other hand, stars that 
have lost all, or most, of the H-envelope produce CC-SNe that are known as "stripped- 
envelope" (briefiy "stripped") SNe. These include, in a sequence of increasing envelope 
stripping. Type lib (He rich, but still showing some H), Type lb (He- rich, no H), and 
Type Ic (deprived of both H and He). Some SNeIc (hereafter "broad-hned" SNeIc) show 
very broad absorption features in optical spectra obtained within a few weeks of the 
explosion; these features are produced by material moving at f > 0.1c (c is the speed of 
light), probably the result of an explosion characterized by a kinetic energy {E^) larger 
than the canonical value of ~ 10^^ erg (3). The most energetic broad-lined SNeIc can 
reach > 10^^ erg (hereafter "hypernovae" or "GRB-HNe") (4), and can be associated 
with gamma-ray bursts (GRBs) (5). Figure 1 summarizes the relation between E^ and 
the mass of ^^Ni that powers the optical light of stripped CC-SNe (6). 

An important unsolved question concerns how gravitational energy of the collapse 
is turned into outward motion of the SN explosion. Most recently proposed scenarios 
involve aspherical explosions {7-11). Therefore, mapping the explosion geometry can 
be illuminating. It is especially critical to establish whether the explosion geometry is 
similar for all CC-SNe, or at least for different subclasses (GRB-HNe, broad-lined SNeIc, 
normal-energy SNe, or stripped vs. Type II SNe). 

Some evidence for asymmetric explosions was obtained from the polarization detected 
in several SNe II {12) and SNe Ib/c {13), and in the broad-lined SN Ic 2002ap {I4-I6). 
However, as there are still few such detections, no systematic study exists. 
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The best way to investigate a supernova's inner ejecta geometry is through late-time 
spectroscopy. At t > 200 days after the explosion, expansion makes the density of the 
ejecta so low that optical photons produced anywhere in the ejecta escape without in- 
teracting with the gas. At these epochs, the SN spectrum is nebular, showing emission 
lines mostly of forbidden transitions. Because the expansion velocity is proportional to 
radius of any point in the ejecta, the Doppler shift indicates where the photon was emit- 
ted: a photon emitted from the near/far side of the ejecta is detected at a shorter/longer 
(blueshifted/redshifted) wavelength. The late-time nebular emission profiles thus probe 
the distribution of the emitting gas within the SN ejecta. This strategy is particularly 
effective for stripped CC-SNe, since we can look directly into the oxygen core. 

Analysis of the late-time spectra of the GRB-HN SNIc 1998bw {lit) and of the broad- 
lined SN Ic 2003jd {18) provided evidence that these objects shared a similar, bipolar 
explosion. However, we viewed SN 1998bw on-axis, and SN2003jd sideways. This seems 
consistent with the fact that SN 1998bw was associated with a GRB, while SN2003jd was 
not, and suggests that if SN 2003jd also produced a GRB, this was missed because of its 
orientation {{18), but see {19) for concerns). 

We obtained late-time spectra of stripped CC-SNe to study their morphology and 
quantify their properties. Our data were obtained mostly with FOCAS {20) on the 8.2 m 
Subaru telescope and with F0RS2 on the ESO Very Large Telescope (VLT). Additional 
data are from {21). The strongest emission hue in stripped CC-SNe is [O l] AA6300, 6363. 
Despite being a doublet, it behaves like a single transition if the lines are sufficiently broad 

0.01c), because the red component is weaker than the blue one by a factor of 3 (see 
supporting online text). 

Since at epochs larger than 200 days after the explosion the ejecta are transparent 
to line emission, and radiation transfer is unimportant (see supporting online text), we 
selected spectra obtained at least 200 days after discovery to build an unbiased data set. 
Additionally, we did not include hypernovae discovered through an associated GRB, to 
avoid bias in the viewing angle. The selection procedure and possible biases are discussed 
in {22). Our sample (Tab. SI) is the largest published data set to date of stripped CC- 
SNe at such late epochs. Figure 2 shows the spectra of the 18 SNe in our sample. Among 
them, 13 are presented here for the first time. 

The observed [O l] AA6300, 6363 emission profiles can be compared with the prediction 
of various explosion models. We use three representative models from {23): one extremely 
aspherical (BP8), one mildly aspherical (BP2), and one spherical (BPl). 

In the spherical model, ^^Ni is confined in a central high-density region with an inner 
hole, and is surrounded by a low-density, oxygen- rich region {2^. This results in a single- 
peaked, but flat-topped [O l] proflle, independent of the orientation. 

On the other hand, the bipolar model {2^-26) is characterized by a low-density ^^Ni- 
rich region located near the jet axis, where the jets convert stellar material (mostly oxygen) 
into Fe-peak elements, and by a high-density disk-like structure composed of unburned 
oxygen- rich material, as the jet expands laterally only weakly {17, 24) (see Fig. SI). The 
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[O l] profile in a bipolar model depends on both the degree of asphericity and the viewing 
angle {23, 27) (sec supporting online text and Fig. SI). If a bipolar supernova explosion 
is viewed from a direction close to the jet axis, the 0-rich material in the equatorial region 
expands in a direction perpendicular to the line of sight, and the [Ol] emission profile is 
observed to be sharp and single-peaked. On the other hand, for a near-equatorial view 
the profile is broader and double-peaked. The best fit to the light curve and the spectra 
of the GRB-HN SN 1998bw was obtained with model BPS {23, 27). 

For the degree of asphericity of this model, the [O l] profile should switch from single- 
to double-peaked near a viewing angle 9 fa 50° measured from the jet direction (Fig. S2). 
The predicted frequency of double-peaked [O l] is thus ~ 64%, in the absence of bias in the 
orientation. For a less aspherical model the fraction of double peaks is reduced: model 
BP2 shows double peaks only for 6 > 70°, and has a double-peak fraction of ~ 34%. 
With this variety of [Ol] profiles, statistics (Tab. 1) can be used to constrain the degree 
of asphericity and remove the uncertainty in the viewing angle. 

Figure 3 shows the [O l] emission profiles for our sample. Out of 18 SNe, 5 (SNe 2003jd, 
2004fe, 2005aj, 2005kl, and 2006T) clearly show double-peaked [Ol] profiles. Four others 
(SNel997dq, 2004gn, 2005kz, and 2005nb) are apparently transitional objects with flat- 
topped or mildly peaked [O l] profiles, and there is a marginal detection of double peaks in 
some cases. The remaining 9 SNe exhibit single-peaked profiles. For illustrative purposes, 
the observed profiles are compared with model predictions in Figure 3. 

The line profiles are well reproduced by the bipolar explosion model assuming dif- 
ferent viewing angles. Although detailed fits are not unique for individual objects, this 
uncertainty does not affect the bulk statistics (single- vs. double-peaked): the presence of 
double-peaked [O l] profiles is not predicted in spherical models, and their fraction yields 
a secure estimate of the number of aspherical SNe viewed sideways, assuming that the 
sample is unbiased in orientation. The high incidence of double-peaked profiles is an im- 
portant discovery: double-peaked [Ol] was previously reported only for SNIc 2003jd {18) 
and SNlb 2004ao {28). The observed fraction of double-peaked profiles, (5 to 9)/18 (= 
28%-50%, median ~ 39%) is consistent with all stripped CC-SNe being mildly aspherical, 
hke model BP2. 

On the other hand, the observed fraction does not support the possibility that all 
stripped CC-SNe are extremely aspherical explosions like SN 1998bw (i.e., model BP8), 
as then we would expect an even larger fraction of double-peaked [O l] profiles: according 
to our Monte Carlo simulations with randomly oriented viewing directions, the observed 
number should be 7 — 16 out of 18 SNe with 99% confidence level and 10 — 13 with 70% 
confidence level. Alternatively, about half of all stripped CC-SNe may have asphericity 
as large as that of GRB-HNe, with a double-peak incidence of ~ 64% (BP8), with the 
remaining half being approximately spherical and yielding the bulk of the single-peaked 
profiles. With a larger SN sample we could look in more detail at [O l] profiles as functions 
of the degree of asphericity and orientation, and more fully explore these scenarios. 

Although our sample is still small, we can look for statistical differences between 
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GRB-HNe, broad-lined SNeIc, and other stripped CC-SNe. In our sample, six SNe are 
broad-lined (non-GRB) SNeIc (SNe 1997ef, 1997dq, 2002ap, 2003jd, 2005kz, and 2005nb), 
while the rest are probably "normal" stripped CC-SNe. The observed fraction of double- 
peaked [Ol] is 33% (1 out of 3 SNe) for broad-hned SNelc and 36% (4 out of 11) for 
the others, if SNe with transitional [Ol] profiles (Fig. 3) are excluded. Within the 
statistical uncertainty, caused especially by the small sample of broad-lined SNe Ic, there 
is no difference between the two groups. Both arc consistent with the predictions of 
model BP2, and show too few double-peaked profiles for model BPS. This suggests 
that on average, broad-lined, non-GRB SNeIc are less aspherical than GRB-HNe, and 
morphologically closer to normal stripped CC-SNe. 

Broad-hned, non-GRB SNeIc have typically smaller £'k than GRB-HNe (Fig. 1). 
There has been speculation that broad-lined, non-GRB SNeIc might be intrinsically sim- 
ilar to GRB-HNe but viewed off-axis, leading to apparently smaller i?K- Although this 
may still be true for a small subset of them, the moderate asphericity inferred for this 
group suggests that broad-hned, non-GRB SNeIc are probably intrinsically different from 
GRB-HNe. 

On the other hand, asphericity is not a special feature of GRB-HNe, but rather a 
generic property of stripped CC-SNe. Both broad-lined and normal stripped CC-SNc 
have a moderate degree of asphericity. All stripped CC-SNe probably share to some 
extent a common explosion mechanism that generates the same kind of asphericity, with 
GRB-HNe likely the most aspherical. Our result offers an important insight into the 
theory of SN explosions. The most popular models for high- and low-energy CC-SNe 
are in fact different: black hole formation and the production of a jet in hypernovae and 
perhaps broad-lined SNeIc {29); delayed neutrino heating from the proto-neutron star for 
normal SNe (30). In the former case, the explosion is probably initiated along the axis of 
rotation and/or magnetic field {10, 11, 29), while in the latter case some asphericity may 
be generated by hydrodynamic instabilities {7-9). The result supports recent theoretical 
scenarios of the supernova explosion, which suggest that an important role in the collapse 
is played by hydrodynamic instability, rotation, or magnetic fields. 

Recently, Modjaz et al. {31) showed another sample of late-time spectra of stripped 
CC-SNe. Their conclusions are similar to ours. One SN in their sample (SN 2004ao with 
double-peaked [Ol]) could be added to our sample according to our selection criteria. It 
increases the frequency of double-peaked events, but it does not change our conclusions 
within the uncertainties. 
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Figure 1: Relation between the kinetic energy of the explosion (-Er) and the mass of 
ejected ^^Ni [M(^^Ni)] of stripped CC-SNe (see supporting online text). Colors indicate 
the ejecta mass (Mej). For SN lib 1993J and SN lb 2005bf, the ejecta mass after sub- 
tracting the He envelope mass is shown as Mgj, to compare with SNe Ic, which lack the 
He envelope. SNe associated with GRBs or an X-ray flash (XRF) are indicated by stars, 
and broad-lined SNe Ic without GRBs/XRFs by open circles. Other (normal) stripped 
CC-SNe are shown as dots. 
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Figure 2: Nebular spectra of the supernovae in the sample used in this study. Narrow 
lines (e.g., Ho; at 6563 A) originate from a diffuse superposed H II region, not from a 
SN. The spectra of SNe 1991N, 1997dq, 1997ef are from (21). The other spectra were 
obtained with the Subaru telescope, except for SN 2006F which was taken by the VLT 
{22). Spectra were deredshifted using the redshift obtained from the observed wavelength 
of the narrow Ha emission if possible; otherwise, the redshift of the nucleus of the host 
galaxy was adopted. For presentation, the flux is arbitrarily scaled and shifted vertically. 
The strongest emission line (dashed line) is [Ol] AA6300, 6363. The feature at ~ 7, 300 A 
is [Ca ll] AA7291, 7324 contaminated by several emission lines: [Oil], [Fe ll], and [Ni ll]. 
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Figure 3: Observed [Ol] AA6300, 6363 emission- line profiles (black curves), classified into 
characteristic profiles: single peaked (denoted by "S," top 9 SNe), transition ("T," middle 
4 SNe), and double peaked ("D," bottom 5 SNe). For presentation, model predictions of 
the bipolar model (23) with different viewing directions are shown for Model BPS (red 
curves, with the direction denoted by the red text) and for the less aspherical model BP2 
(blue). The models shown here have smaller expansion velocities (corresponding to E^^ 
a few xlO^^ erg) than the one applied to SN 1998bw {Ek > 10^^ erg) (23). The separation 
between the two peaks in the double-peaked cases is ~ 0.01c — 0.02c, much larger than 
the wavelength resolution in the observations (~ 10 A, corresponding to 0.0015c). 
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Table 1: Fraction of double-peaked [Ol] SNe. 



Model 


Dividing angle^ 


Fraction^ 


Spherical 




0% 


BP2 


- 70° 


~34% 


BPS 


~ 50° 


- 64% 


Observed 




39 ± 11% 



^The viewing angle (^o measured from the jet direction) which divides the [Ol] profile — 
i.e., single peaked if < and double peaked if ^ ^ ^ 90°. 
^Fraction of SNe showing double-peaked [Ol]. 
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SUPPORTING ONLINE MATERIAL 



Asphericity in Supernova Explosions from Late-Time Spectroscopy 

Keiichi Maeda et al. 

Data and Methods 

Sample of Supernova Nebular Spectra: 

In the years 2005 to early 2007, we obtained a number of spectra of stripped CC-SNe using 
the 8.2 m Subaru telescope equipped with FOCAS. We also obtained several spectra using 
the VLT with F0RS2 from late 2006 to early 2007. Some of the spectra are distinctly 
nebular, while others are not. In addition to this new sample, we include some spectra 
presented in (SI), which contains the largest sample of previously published late-time 
spectra of stripped CC-SNe. The spectral resolution is ~ 10 A for the spectra obtained 
by the Subaru or VLT. The spectra from the previous sample (SI) were obtained with 
the Shane 3-m reflector at Lick Observatory, with similar spectral resolution (6-15 A). 

Because of the faintness of SNe in the late phases, the spectra can be obtained only for 
relatively bright SNe even with the largest telescopes in the world. The limiting functions 
are thus the intrinsic luminosity and the distance. The luminosity bias may be important 
when interpreting the data if the effect is strong: our sample then represents a subclass 
of luminous SNe (i.e., SNe with a large amount of ^^Ni), and may affect the asymmetry 
statistics we investigate in the present paper. 

However, we believe that the effect is not strong enough to affect our interpretation. 
First, our criterion for successful spectroscopy is ~ 23 mag. This roughly corresponds to 
SNe with a peak magnitude of 17 or 18. This brightness is similar to the criterion for 
obtaining high-quality data with typical SN observations in the early phases using smaller 
telescopes. As such, our sample should not be very different from samples of well-studied 
stripped CC-SNe in past studies using early-phase observations which, for example, led to 
the result shown in Figure 1. Second, the peak luminosity of very luminous stripped CC- 
SNe is larger than that of typical ones by at most a factor of ~ 5 (e.g., SN 1998bw {S2)), 
so the detection volume of these luminous events is larger than that of typical ones by at 
most one order of magnitude. However, the intrinsic rate of the very luminous events is 
probably smaller than that of the typical ones by nearly two orders of magnitiide {S3). 
Thus, contamination by very luminous SNe should not be large in our sample. Finally, 
there is a correlation between the ^^Ni mass and the early-phase width of absorption lines 
{S4)- Using this relation, we indeed divide our sample into a potentially luminous class 
(i.e., broad-lined SNe Ic) and the other ones in the discussion of the asphericity statistics. 
Also, according to past studies, normal stripped CC-SNe (i.e., not broad-lined SNe Ic) 
have similar ^^Ni mass (~ O.IMq) and represent a wide range of progenitor masses up 
to ~ 2OM0 {S4)- This suggests that the luminosity bias, if present, is not strong in the 
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normal stripped CC-SNe. In summary, the luminosity bias should not be strong in our 
sample of late-time spectra. Thus, the distance is probably the limiting factor, and our 
data consist of SNe within a distance of at most ~ 100 Mpc. 

Other possible biases are (1) a possibility of preferentially observing SNe from a special 
direction (e.g., a jet axis), and (2) contamination of SNe which are young and dense so 
that their ejecta are not sufficiently transparent to trace the geometry by the [O l] profile. 
Here we describe how we selected the final unbiased list of supernova spectra expected to 
be at a fully nebular epoch to study the intrinsic line profile (Table SI). 

First, we excluded SNe 1998bw, 2003dh, 2005bf, and 2006aj. SN Ic 2006aj was associ- 
ated with an X-Ray flash (XRF: a low-energy analog of a GRB) {S5); thus, the detection 
of SN 2006aj is probably biased by the viewing-angle effect, being easier for an observer 
along the jet direction. For the same reason, we excluded SNe 1998bw and 2003dh, associ- 
ated with GRB 980425 and GRB 030329, respectively, although there are some published 
late-phase observations for these objects {S6, ST). SN lb 2005bf is reported to be a very 
peculiar SN, whose explosion mechanism may be different from that of typical CC-SNe; 
the nebular observations are found to be peculiar, too {S8). Observations and analyses 
of SN 2006aj and SN 2005bf in nebular phases arc described in separate papers {S8-10). 
Note that in the late-time phases considered here, the dependence of luminosity on the 
viewing angle is negligible. 

Next, we excluded supernovae which were observed only at times before entering the 
fully nebular phase, since optical transport effects can change the hue profile. The choice 
is not trivial; the optical depth (r) differs for different supernovae (r oc M^/£'k), and 
the ejecta mass (Afcj) and the ejecta kinetic energy (-Ek) are not known a priori. The 
uncertainty in the explosion date further complicates matters. 

Therefore, we apphed a very simple criterion. We checked the phases of the obser- 
vations with respect to the discovery date (t), and selected a supernova only if it was 
observed at t > 200 days at least once. Thus, any SNe possibly younger than 200 days 
were excluded from the list. This criterion, t > 200 days, is sufficiently late so that any 
selected supernova is certainly in the nebular phase. 

Following these procedures, we obtained the unbiased sample of stripped CC-SNe 
whose spectra are distinctly nebular. The sample contains 18 SNe. One of the SNe is of 
Type lib. The sample contains 13 SNe newly observed (11 only by the Subaru telescope, 
and 2 both by the Subaru and VLT). Spectra of two broad-lined SNe Ic (2002ap and 
2003jd) were also obtained with the Subaru telescope, which were already reported {Sll, 
SI 2). The remaining 3 SNe are from the previously published sample of {SI). In terms 
of the spectral characteristics in the early phase, 6 SNe in our hst have been classified 
as broad-lined SNe Ib/c without association with a GRB (SNe 1997dq, 1997ef, 2002ap, 
2003jd, 2005kz, and 2005nb - see also Figure 1 in the main text). 

In cases where the supernova was observed at multiple epochs, we included the object 
in the nebular sample if the last observation satisfied the condition t > 200 days. Such 
SNe were used to check the vahdity of the selection criterion. SN 2005aj was observed 
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at three epochs {t — 189, 249, and 312 days), showing significant evolution of the [Ol] 
AA6300, 6363 profile between the first two epochs. SN 2006F was observed at t = 173 and 
312 days, showing no noticeable evolution in the [Ol] profile. The other SNe observed 
multiple times at t > 200 days do not show noticeable evolution in the [O l] profile. These 
observations justify our criterion {t > 200 days). 

All the spectra shown in Figures 2 and 3 in the main text are at t > 200 days except 
for that of SN 2005nb. The signal-to-noise ratio for the spectrum of SN 2005nb at t — 403 
days is low, but the [O l] profile does not evolve significantly from t — 195 days (shown in 
Figures 2 and 3 in the main text). 

Supporting Text 

The properties of stripped CC-SNe: E^, M(5^Ni), and Mej: 
The kinetic energy of the explosion (Ek), the mass of ^^Ni [M(^^Ni)], and the ejecta 
mass (Mej) shown in Figure 1 of the main text are derived by modeling the early-phase 
observations using approximate analytical expressions (SNe 2003jd and 2004aw) or one- 
dimensional radiation transport codes (the other SNe). An exception is SN lb 2005bf, for 
which M(^^Ni) is derived by modeling the late-time observations {S8). 

Assuming spherical symmetry, the estimate of Ek (shown in Figure 1) represents an 
isotropic value, which could be different from the intrinsic value by up to a factor of a few. 
For example, for SN 1998bw (associated with GRB 980425), we estimated Ek ^ 5 x 10^^ 
erg with a spherical model {S2), but Ek ?a 2 x 10^^ erg with a two-dimensional jet model 
{SI 3). The difference can be attributed to a viewing-angle effect. Since SN 1998bw is an 
extreme case (see the main text), the difference between the isotropic value and the real 
one should be smaller for other SNe. 

The references are as follows: SN lib 1993J {SI4), SN Ic 19941 {S15), SN Ic 1997ef 
(SI 6), SN Ic 1998bw {S2), SN Ic 2002ap (5i7), SN Ic 2003dh {SI 8), SN Ic 2003jd {SI 9), 
SN Ic 20031w {S20), SN Ic 2004aw {S19), SN lb 2005bf {S8, S21), and SN Ic 2006aj {S22). 

The Aspherical SN Model: 

The aspherical models presented in this paper are from {S23, S24)- The jet-driven explo- 
sion is simulated by promptly depositing the kinetic and thermal energy in the innermost 
core. The initial kinetic energy is distributed in a jet-like way, applying Vz = olz and 
Vr = /3r, where (r, z) is a cylindrical coordinate system (assuming axial symmetry with 
respect to the z-axis), and and Vr are the initial deposited velocity components. Mod- 
els BP8 and BP2 assume a/P = 8 and 2, respectively; see {S23, S24) for details. The 
distributions of density, ^^Ni, and ^^O are shown in Figure SI. 

The optical depth to the [O l] A6300 can be approximated by 0.01Mo/Mo(tcxp/200 days)~^ 
for a homogeneous medium with an oxygen mass of Mq and an expansion velocity of 0.02c 
(where texp is time since the explosion) {S25). The [O l] A6363 line is always more trans- 
parent than the 6300 A component because the latter has a larger transition probability. 
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Thus, the effect of the optical depth to the [O l] hne profile should be neghgible. Fur- 
thermore, wc have confirmed that the [O l] line is always transparent in om^ models at 
tcxp > 200 days. In this optically thin case, the blue (6300 A) component is stronger than 
the red (6363 A) one by a factor of 3, and the doublet behaves like a single transition. 

The jet-driven model predicts unique features in the [O l] profile (Figs. SI and S2) 
{S23, S24)- The SN ejecta kinematics are described by homologous expansion, R — 
v{R)tcxp, where R and v{R) arc respectively the radius (in a spherical coordinate) and 
the radial expansion velocity of a point in the ejecta. Non-radial velocity components 
should be negligible. The observed wavelength of a photon depends on the point from 
which it was emitted, as A = Ao(l — v\\/c), where Aq is the rest wavelength of the line, 
is the line-of-sight velocity toward the observer, and c is the speed of light. Then, defining 
d as the projection of it! onto the line of sight, it can be shown that = d/tcxp, and this 
is the same in a given plane perpendicular to the line of sight. This means that, in the 
late phase, all photons emitted at the same depth along the line-of-sight have the same 
wavelength. The wavelength of a photon emitted from the near/far side of the ejecta is 
detected as shorter /longer (blueshifted/redshifted). Thanks to these characteristics, and 
to the transparency of the ejecta, the late-time emission profile can be used to probe 
the density distribution of the emitting gas within the SN ejecta (the "SN scan"), just 
like Computerized Tomography scans the human body (the "CT scan"). The SN scan is 
shown in Figure SI for observers along the z-direction and r-direction, respectively. The 
model prediction (BPS) for the [O l] profile is shown in Figures SI and S2 as a function 
of the viewing angle {9). 

If we view the bipolar supernova explosion from a direction close to the jet axis, most 
of the 0-rich material is in the equatorial region, expanding in a direction perpendicular 
to the line of sight and producing a narrow-peaked [Ol] emission profile (Fig. S2). The 
SN scan provides a monotonically increasing amount of 0-rich material as one probes 
nearer the center of the SN, then a monotonically decreasing amount of 0-rich material 
as one moves away from the center (Fig. SI). On the other hand, the [Ol] line shows a 
double-peaked profile if viewed from near the equatorial direction. In this case, the SN 
scan results in an increasing amount of 0-rich material as one moves from the observer 
to the edge of the inner hole of the 0-rich disk, then a decrease toward the center of the 
SN, where the non-emitting hole occupies the largest area. The amount of O then again 
increases as one moves from the center to the edge of the hole on the other side, and 
subsequently it again decreases toward the receding side of the ejecta (Fig. SI). Thus, 
two characteristic Doppler-shifted emission peaks (blueshifted and redshifted) appear. 
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|(b. d) 




Figure SI: The bipolar Model BPS {S23, S24). (a) Distribution of ^^Ni (which decays 
into cobalt and then iron, shown in blue) and oxygen (red). Density contours are shown 
covering two orders of magnitude and divided into 10 equal intervals on a logarithmic 
scale, (b) The "SN scan" for the [O l] line as seen from an observer placed near the axis 
of the jet (z). (c) The SN scan for [Ol] for an observer in the equatorial direction {x — y 
plane), (d) The theoretical [O l] profile for an observer along the z-direction, and (e) along 
the equatorial direction. 
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Figure S2: Predicted [Ol] line profiles in the aspherical, jet-driven models BPS (red) and 
BP2 (blue) as a function of the viewing angle 9 {S24)- 
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Table SI: Supernova Samples 



SN 


Type 


Observing Date 


Days"*^ 


References 


Figure^ 


1991N 


Ic 


1992-1-9 


286 


Lick {SI) 


o 


1997dq 


broad-Ic 


1998-8-18 


289 


Lick (SI) 


o 


1997ef 


broad-Ic 


1998-9-21 


299 


Lick (SI) 


o 


2002 ap 


l)roa(l-Ic 


2002-9-15 


229 


Sul)aru (5i2) 


o 


2003jd 


broad-Ic 


2004-9-12 


323 


Subaru {Sll) 


o 


2004dk 


lb 


2005-7-6 


342 


Subaru 








2005-8-26 


392 


Subaru 


o 


2004fe 


Ic 


2005-7-6 


250 


Subaru 


o 






2005-8-26 


300 


Subaru 




2004gk 


Ic 


2005-7-6 


224 


Subaru 


o 


2004gn 


Ic 


2005-7-6 


218 


Subaru 


o 


2004gq 


lb 


2005-8-26 


258 


Subaru 








2005-10-25 


318 


Subaru 


o 






2005-12-27 


381 


Subaru 




2004gv 


Ib/c 


2005-7-6 


206 


Subaru 








2005-8-26 


256 


Subaru 


o 


2005aj 


Ic 


2005-8-26 


189 


Subaru 








2005-10-25 


249 


Subaru 








2005-12-27 


312 


Subaru 


o 


2005kl 


Ic 


2006-6-30 


220 


Subaru 


o 






2006-12-25 


398 


Subaru 




2005kz 


broad- Ic 


2006-6-30 


211 


Subaru 


o 


2005nb 


broad-Ic 


2006-6-30 


195 


Subaru 


o 






2007-1-24 


403 


Subaru 




2006F 


lb 


2006-6-30 


173 


Subaru 








2006-11-16 


312 


VLT 


o 


2006T 


lib 


2006-11-26 


301 


Subaru 








2006-12-25 


329 


Subaru 


o 






2007-2-18 


384 


VLT 




2006ck 


Ic 


2007-1-24 


249 


Subaru 


o 



^Days since the discovery. 

^The spectra shown in Figures 1 and 2 in the main text are marked by circles. 
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